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Abstract 

We study the resonant x-ray scattering at Si and Al K-edges from chiral materials, a-quartz 
and a-berlinite. We derive the general form of the scattering matrix for the dipole transition by 
summing up the local scattering matrices which satisfy the symmetry requirement. The oscillation 
term is obtained in the spectral intensity as a function of azimuthal angle with an expression of 
possible phase shift. We evaluate the parameters undetermined by the symmetry argument alone 
on the basis of underlying electronic structures given by the bond-orbital model. The spectra are 
calculated on forbidden spots (001), (OOl), (002), and (002) in circular polarizations without ad- 
justable parameter, reproducing well the experimental curves depending on polarization, chirality, 
and scattering vector. Some discrepancies remain in the phase shift in a-quartz. 

PACS numbers: 61.05.cc, 71.20.Nr, 78.70.Ck 
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Sample 

FIG. 1: Scattering geometry. The wave vectors kj and kj are for the incident and scattered photons, 
respectively. The scattering vector is defined by G = kj — kj. The e^, g-k are polarization vectors 
for the incident photon, and and are for the scattered photon. The Cq-, ej^ are perpendicular to 
the scattering plane. Three vectors (eg-, e,ri kj) and (ej^, e^, kj) constitute a right-handed coordinate 
frame, respectively. The sample is right-handedly rotated around G by azimuthal angle ^. 

I. INTRODUCTION 

The a-quartz (Si02) and a-berlinite (AIPO4) are known to have two crystal forms, the 
right-handed screw (space group No. 152, P3i21) and the left-handed screw (No. 154, P3221). 
One is the mirror image of the other, thus the crystals are called to have different chirality. 
The two forms have been distinguished by using the optical activity since the discovery 
of Arago and Biot;- as linearly polarized light passes through the crystal, the direction 
of polarization rotates about the beam axis oppositely according to the chirality. Another 
method to distinguish chirality is the anomalous x-ray scattering by which atomic scattering 
amplitudes become complex numbers,-*^ leading to the determination of atomic positions for 
systems with different chirality. Note that the conventional x-ray diffraction using Thomson 
scattering could not distinguish chirality. 

Recently, another method has been attempted to distinguish chirality: the resonant x- 
ray scattering (RXS) with circularly polarized beam.->^ The RXS process may be described 
at the Si K-edge in a-quartz and at the Al K-edge in a-berlinite as follows. The Is-core 
electron is prompted to unoccupied p-symmetric states by absorbing photon [electric-dipole 
(El) transition], and subsequently the excited electron is recombined to the core-hole by 
emitting photon (El transition). This will be called the El-El process. Figure [T] shows 
the scattering geometry, where the sample is rotated by azimuthal angle \I/ around the 



scattering vector G = kj — kj. It is known that the tensor character of the scattering matrix 
of RXS could give rise to the intensity on the spots forbidden in Thomson scattering.-""- 
Measuring the spectra on the forbidden spots G = (001) and (OOl), Tanaka et alS^ have 
found characteristic patterns depending on chirahty in the spectral intensity by means of 
switching polarizations from the right-handed one (RCP) to the left-handed one (LCP). 

The purpose of this paper is to analyze such spectra from underlying electronic structures. 
We start by deriving the scattering matrix on the basis of symmetry requirement; the only 
assumption made is that the total resonant scattering matrix is the sum of the local scattering 
matrix" on each Si or Al site. We introduce the local dipole-dipole correlation function 
instead of the local scattering matrix in this procedure, since the former quantity consisting 
of real numbers makes the expression transparent. We obtain a general formula of scattering 
matrix depending on polarization, chirality, and scattering vector. The scattering intensity 
contains the oscillation term as a function of azimuthal angle, which has a phase shift and 
the amplitude with chirality dependence different from previous ones.-'^'"— 

The general expression of scattering matrix involves parameters undetermined by the 
symmetry argument alone. In this paper, we evaluate these parameters from underlying 
electronic structures by exploiting a simple bond-orbital model developed by Harrison.— 
This model is known to work well on the ground-state properties in covalent-bonding sys- 
tems; it considers the strong coupling between the sp^-hybrid on Si atoms and the sp^'^^- 
hybrid on O atoms with four bonds per Si atom in a-quartz. In a-berlinite. Si atoms are 
replaced by Al and P atoms. In the ground state, the bonding states are occupied, while 
the anti-bonding states are unoccupied. In the intermediate state, one anti-bonding state 
is occupied at the core-hole site, on which the core-hole potential works. The parameters 
in the scattering matrix are evaluated by using the electronic structures thus determined. 
The spectral intensities are calculated depending on polarization, scattering vectors, and 
chirality as a function of azimuthal angle, reproducing the experimental curves particularly 
well for a-berlinite. The present model calculation predicts the phase shift of oscillation to 
be TT for both a-quartz and a-berlinite. The experimental curves indicate the phase shift vr 
in a-berlinite and its deviation from tt in a-quartz.^'^ Since both materials are expected to 
have similar electronic structures, such a difference is puzzling to us. The possible origin for 
the deviation will be discussed in the last section. 

The present paper is organized as follows. In Sec. HIl the crystal structures of a-quartz 
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and a-berlinite are briefly described. In Sec. Illlt the scattering matrix is formulated from 
the symmetry requirement. In Sec. IIVI the RXS intensity as a function of azimuthal angle 
is formulated for the incident x-ray beam specified by the Stokes parameters. In Sec. El the 
bond-orbital model is introduced to evaluate the electronic structures as well as the RXS 
intensity. In Sec. Wl\ the calculated results are discussed in comparison with experiments. 
Section IVIII is devoted to the concluding remarks. 

II. CRYSTAL STRUCTURE 

A. a- quartz 

The crystal structure of a-quartz is hexagonal with three Si atoms per unit cell with 
a = b = 4.91A, and c = 5.40A. It has chirality described by two different space group P3i21 
(No. 152), the right-handed screw, and P3221 (No. 154), the left-handed screw. As shown in 
Fig. Si atoms sit at the positions {u, 0, 0), (0, u, |), (1 — m, 1 — u, |) for No. 152, while at 
the positions (u, 0, 0), (1 — m, 1 — u, |), (0, u, |) for No. 154, where u = 0.47. One O atom sits 
between each nearest-neighboring Si-Si pair. Each Si atom is surrounded by a tetrahedron 
of O atoms. We introduce the Cartesian frame where x and z axes are parallel to the 
crystal a and c axes and the origin is set at the center of the Si atom labeled as 1 in Fig. 
[2J Then, the coordinates of O atoms of the tetrahedron are given by (0.93, =fO. 59, —1.13), 
(-0.96, 1.13, TOM), (-0.96, -1.13, ±0.64), (0.93, ±0.59, 1.13) in units of A, where the upper 
and lower signs correspond to No. 152 and No. 154, respectively.^^ Note that there exists the 
symmetry of two- fold rotation around the a-axis for both No. 152 and No. 154, and that 
the crystal of No. 154 is the mirror image of No. 152 with respect to the a-b plane. The 
settings correspond to the r (±) and z (— ) settings for the crystals of No. 152 and No. 154, 
respectively, in which the twofold axis develops (±) charge at its positive end on crystal 
extension along the axis. Other tetrahedrons are given by rotating ^ (~^) ^ind ^ (~^) 
around the c axis with translating |c, |c along the c axis for No. 152 (No. 154). 

B. a-berlinite 

The crystal structure is close to the a-quartz structure, where Si atoms are replaced by 
Al and P atoms alternatively along the c-axis. Therefore, the unit cell is doubled along 
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(a) No. 152 (b) No. 154 




FIG. 2: Si atoms projected onto the a-b plane in the unit ceh of a-quartz. Atoms labeled by 1, 2, 
3 are located at {u, 0, 0), (0, u, |c), {1 — u,l — u, |c) for No. 152, and at {u, 0, 0), (1 — u, 1 — u, |c), 
(0, u, |c), respectively, u = 0.47. Si atoms are surromided by O atoms forming a tetrahedron. 

the c-axis; there are three Al atoms and three P atoms per unit cell with a = b = 4.94A, 
c = IO.95A. More precisely, for No. 152, Al atoms sit at (u, 0, 0), (0, u, |), (1 — m, 1 — u, |), 
and P atoms at (1 -u',1 - u', |), {u', 0, |), (0, u', |), with u ^ u' = 0.47, while for No. 154, 
Al atoms sit at (^,0,0), (1 — m, 1 — m, |), (0,u, |), and P atoms at (0,u', |), (m',0,|), 
(1 — n', 1 — m', |). In the Cartesian frame where the x and z axes are along the a and 
c axes, respectively, and its origin is the center of the Al site at (m,0,0), the coordinates 
of the O atoms of the tetrahedron are given by (1.01, ±0.67, —1.25), (—0.95, 1.27, ±0.69), 
(—0.95, —1.27, =f0.69), (1.01, =f0.67, 1.25) in units of A, where the upper (lower) sign corre- 
sponds to No. 152 (No. 154).— Note that the tetrahedron surrounding Al for No. 154 is similar 
to that for No. 152 in a-quartz. Other tetrahedrons surrounding Al are given by rotating 
T T (^t) s-^ou^d the c-axis with translating |c, |c along the c-axis for No. 152 

(No. 154), according to the screw symmetry. 



III. SCATTERING MATRIX 



Let the incident and scattered photon polarizations be specified as and Xa in Cartesian 
frame where the xi(= x) and a;3(= z) axes are along the a and c axes. The scattering matrix 
at site j, M{j; u), may be expressed as 

[MO;.)U^ <f»<"iy^ . (3.1) 

where the dipole operator Xa is measured from the center of site j. Ket \g) represents 
the ground state with energy eg, and \n) represents the intermediate state with energy e„ 
(including the core- hole energy). The F represents the life-time broadening width of the 
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Is- core hole. 

The total resonant scattering matrix could be well approximated by the sum of the local 
amplitudes at Si sites or at Al sites, since the Is state is localized at each Si site. The 
scattering geometry is shown in Fig. [T] For the incident and scattered wave vectors kj and 
kj, the total resonant scattering matrix may be expressed as 

[M(G,a;)],,^ = ^[M(j;u;)],,^exp(-zG ■ r,). (3.2) 

j 

where j runs over Si sites for Si i^-edge (quartz) and over Al sites for Al i^-edge (berlinite). 
The scattering vector is defined by G = kj — kj. 

To analyze the symmetry of the scattering matrix, it may be convenient to introduce the 
local dipole-dipole correlation function at the core-hole site j, which is defined by 

[p(j;e)]a,/3 = ^{9\xa\n){n\xp\g)5{e + eg - en), (3.3) 

n 

and the total dipole-dipole correlation function associated with G, which is defined by 

p(G; e) = ^ p(j; e) exp(-^G ■ r,). (3.4) 

j 

Using the latter quantity, we may express the total resonant scattering matrix as 

M{G;co)= [ ^^^''\ de. (3.5) 
J u - e + iL 

Let the local dipole-dipole correlation function at site (m, 0, 0) be p'^^ with + and — signs 
corresponding to No. 152 and No. 154. It should take the following matrix form according to 
the local symmetry: 



(3.6) 



b{e) ±d{e) 
y ±die) c{e) J 

where a(e), 6(e), c(e), d{e) are real functions. The presence of the off-diagonal elements is 
due to the lack of the inversion symmetry around the Si atom in a-quartz or the Al atom in 
a-berlinite. The zero components are originated from the symmetry of the two-fold rotation 
around the a-axis. The ± signs are originated from the mirror-image relation with respect 
to the a-b plane between No. 152 and No. 154, respectively. 
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Now we define the local correlation function rotated by ±27r/3 around the c-axis: 

fi(2V3)pf'(e)fi-i(2,r/3), (3.7) 
fl(-2ir/3)(iJi*'(e)fl"'(-2r/3). (3-8) 



where rotation matrix is defined by 



R{e) 



( 



cos 6 
sin 6* 




sine 
cosO 
1 



(3.9) 



For the crystal of No. 152, the local dipole-dipole correlation function at (0,m ^c) and at 



(1 — M, 1 — u, |c) are given by p\^\e) and p^^i{e), respectively. On the other hand, for the 
crystal of No. 154, the local dipole-dipole correlation function at {1 — u,l — u, he) and at 



(0,M, |c) are given by p^_i{e) and p\ '{e), respectively. 



1. G = (001) 



Summing up the local correlation function defined above with weight exp{—iG ■ rj), we 
obtain the total dipole-dipole correlation functions for No. 152 and No. 154 as 

p((001); e)i52 = pi^\e) + p[^\e) exp(-227r/3) + p^^,\e) exp(227r/3) 



3 
4 



^ [a{e)-b{e)] i[a{e)-b{e)] 2id{e)^ 



\ 



i[a{e)-b{e)] -[a(e) -6(e)] 2d{e) 
2id{e) 2d{e) J 

p((001);e)i54 = A~\e)+p[-\e)exp{t27r/3)+th^)^M-^'^^/^) 
I \a{t)-h{t)] ~i\a{t) -h{t)] 2id{t)^ 

3 



(3.10) 
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V 



[a(e) -6(e)] -i[a(e) - 6(e)] 2id{e) 
-i\a{e) - 6(e)] -(a(e) - 6(e)) -2rf(e) • (3.11) 
2id{e) -2d{e) j 

Note that there exist only two independent components a(e) — 6(e) and d{e) in the matrix. 
Hence, from Eq. (13.51) . the total resonant scattering matrix is written as 



M((001);w) 



^ A{io) ±iA{uj) tB{uj)^ 



±iA{uj) -A{uj) ±B{uj) 
iB{u) ±B{uj) 



(3.12) 



with 

. / N 3 /■ a(e) — 6(e) , 

where the upper and lower signs correspond to No. 152 and No. 154, respectively. Note that 
74(0;) and B{uj) are complex numbers because of the presence of T. 



2. G = (001) 



The total dipole-dipole correlation functions for No. 152 and No. 154 are given by 

p((OOT); e)i52 = pi''\e) + pt\e) exp(z27r/3) + pL+i^(e) exp(-z27r/3), 
p((OOT); e)i54 = pV{^) + pV{^) exp(-i27r/3) + (e) exp(+i27r/3). 



(3.15) 
(3.16) 



Since a(e), 6(e), c(e), d{e) in Po^^(e) and p±i^(e) are real functions, we notice that the total 
resonant scattering amplitude is obtained from M((001);a;) by replacing ±i by ^i. The 
result is 



M((001);cc;) 



TiA{oj) -A{u) ±B{u) 
y-tB{co) ±B{uj) 
where the upper and lower signs correspond to No. 152 and No. 154, respectively. 



(3.17) 



3. G = (002) and (002) 

Now that the phase factors exp(±i47r/3) are equivalent to exp(^i27r/3), we immediately 
obtain the following relation: 

M((002);a;) = M((OOT);a;), (3.18) 
M((002);a;) = M((001);a;). (3.19) 
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(a) G = (001) 



(b) G = (001) 



b' 




FIG. 3: Top view of the scattering geometry at the azimuthal angle ^' = 0. (a)G = (001): view 
looking down along the c axis from the top of the axis. (b)G = (001): view looking up along the c 
axis from the bottom of the axis. Vectors a and b are translational vectors along the a and b axes, 
respectively. Vector a* is the reciprocal lattice vector conjugate to a. 

IV. POLARIZATION ANALYSIS WITH ROTATING CRYSTAL 

Specifying the polarization vectors by (e^) and e^r (e^) for incident (scattered) photon, 
we formally write the component of the total resonant scattering matrix as 



and so on. We will evaluate these values with rotating crystal in the scattering geometry 
described below. 

A. G = (001) and (002) 

We rotate the crystal right-handedly around G with azimuthal angle Following the 
experimental setup by Tanaka et al.,^^ we define the origin of ^ such that the scattering 
plane contains the b axis, that is, it is perpendicular to the reciprocal lattice vector a* 
conjugate to the translational vector a along the shown in Fig. [3]^a). Note that the 

rotation of the crystal indicates that the scattering plane is inversely rotated with respect 
to the crystal. 

The polarization vectors, which are represented in the Cartesian frame with x and z axes 






(4.2) 
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along the crystal a and c axes, are given by 

= = (-sinf , -cos^',0), (4.3) 

= (— sin cos ^, sin 6*6 sin 'I', — cos 6*3), (4.4) 

= (sin 6'b cos — sin sin 'I', — cos 6*6), (4.5) 

where '5' = — 27r/3, and 9^ is the Bragg angle of the scattering. Using Eq. fl3.12p . we 
obtain 

M((001);u;).v = -Ae^^'*, (4.6) 

M((001);a;)^v = TiA?>me^e^''^^ ± B cosO^e^'^ , (4.7) 

M((001);a;),/^ = ±zA sin ^ee^^^* ± 5 cos ^ee^^*, (4.8) 

M((001);a;)^.^ = -Asin^ 9^0^''^^ , (4.9) 

where the upper (lower) sign corresponds to the crystal of No. 152 (No. 154). We have 
abbreviated A{ijj) and B{ijj) by A and respectively. For G = (002), the corresponding 
components of scattering matrix are given by replacing ±i by =l=z in Eqs. fl4.6p - fl4.9p . 



1. Linear polarization 

For the scattering channels with cr — )■ a' and vr — t- tt', the scattering intensities are the 
same for both No. 152 and No. 154, and are constant as a function of azimuthal angle. On the 
other hand, for the scattering channels with a — )■ vr' and tt — ?■ cr', the scattering intensities 
oscillate function of azimuthal angle. 

For G = (001), the scattering intensities in the o" — t- cr' and a — > tt' channels are given by 

/((001);a;),v = \A\\ (4.10) 
/((001);a;)^v = |Ap sin^ + |5p cos^ T sin 20b sin(3^ ^ 5), (4.11) 

where the upper sign (lower) sign corresponds to the crystal of No. 152 (No. 154), and B*A 
is replaced by lAiJle**^. The phase 5 could take any value in principle. 

For G = (002), the corresponding intensities are given by Eqs. fl4.10p and f l4.1ip with 
replacing sin(3\l' =F 5) by sin(3\t' ± 5). Since 3^ = 3\E' — 27r, we could safely replace by \E'. 
The constant term as well as the amplitudes of the oscillation are the same for both No. 152 
and No. 154. Therefore, the knowledge of the phase shift 5 of the oscillation is necessary in 
order to distinguish the chirality. 
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2. Circular 'polarization of incident beam 



We consider the case where the incident beam is circularly polarized. To include the 
partial polarization, we introduce the polarization density matrix,— which is represented on 
the bases e^^ and e^: 

if 1 + P3 Pi - iP2 \ 
P = 7^\ ], (4.12) 

2 \P,+tP2 I-P3 J 

where Pi, P2, P3 are the Stokes parameters. All three parameters take values between — 1 
and +1. In the unpolarized state. Pi = P2 = P3 = 0; for a completely polarized photon, 

+ P2 + P3^ = 1; P2 = +1(— 1) for the right (left)-handed circular polarization. 

Without analyzing the polarization of scattered x-ray, the scattering intensity may be 
expressed in the following form: 



Pn2ni 1 

n n\ 712 

= f/iPi + U2P2 + Ut\{l + P3) + U^\{1 - P3). (4.13) 

For G = (001), the coefficients are given by 

Ui = ^|AP|cos^B(l + sin2^B)cos(3^T5), (4.14) 

U2 = ±|A|2sin^B(l + sin2^B) + |^5|cos3^Bsin(3^T5), (4.15) 

= |Ap(l + sin2^B) + |5pcos2^BT|^5|sin2^Bsin(3^T5), (4.16) 

= \A\'^sm^9Bil + sm^eB) + \B\'^cos^9B±\AB\sm29Bsm{3^T5). (4.17) 

Thereby the scattering intensity with Pi = is expressed as 

/((OOl); tu) = /o((001); to) + /i((001); to) sin(3^ T S), (4.18) 

with 

/o((001);u;) = |P]^ cos^ ^b + 1^^^^^^ 

X (1 + sin2^B ±2P2sin^B + i'3Cos2^B), (4.19) 

/i((001);a;) = |AP| cos^b(A cos^ ^b T 2P3 sIu^b)- (4.20) 

For G = (002), the corresponding coefficients are given by the same forms of Eqs. f l4.14p - 

( I4.17P with replacing 3\l'=F(5by3\l'±5 and with reversing the signs in Eq. (I4.15p . Hence 
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the scattering intensity with Pi = is expressed as 



/((002); co) = /o((002); co) + /i((002); u) sin(3^ ± 6), (4.21) 

with 

/o((002);u;) = \B]' cos' 6^ + \A\'^-±^^ 

X (1 + sin2 0B T2P2sin0B + ^3Cos2^b), (4.22) 

/i((002);w) = |AP|cos^B(-P2Cos2^BT2P3sin^B). (4.23) 



B. G = (001) and G = (002) 

Following the experimental setup by Tanaka et ai.,-'^ we rotate the crystal by angle vr 
around the reciprocal lattice vector a* in order to use the reverse side of the crystal for the 
scattering, as shown in Fig. |3]^b). The other scattering conditions are kept the same as those 
for G = (001) and (002). Since the direction of G is now opposite to (001), the right-handed 
rotation of crystal around the direction of G means the left-handed rotation around the c 
axis. The polarization vectors, which are represented in the Cartesian frame with the x and 



z axes along the crystal a and c axes, are given by 

= = (sin^l', -cos^f,0), (4.24) 

= (sin 6'b cos ^, sin6'B sin 'i', cos6'b), (4.25) 

= (-sin6'BCOs^', -sin6'Bsini|',cos6'B), (4.26) 

where ^ = + 27r/3. Thereby, from Eq. f l3.17p . we obtain 

M((OOT);a;).v = -Ae^'^'^, (4.27) 

M((OOT);w)^v = TiAsinOBe'^'''^ TBcose^e^'^, (4.28) 

M((OOT);w)^.^ = ii^sin^Be^^^* T^cos^Be^^*, (4.29) 

M((OOT);w)^/^ = -Asin^^Be^^^*. (4.30) 



These equations are nothing but Eqs. ( l4.6p -( lT9ll with reversing signs in front of B. For G = 
(002), M((002); a;)o-'o- ■ • • sue given by the same forms of Eqs. fl4.27p - fl4.30p with replacing 
±i by =Fi. With these expressions, the following results straightforwardly come out. 
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1. Linear polarization 



We obtain the scattering intensity, 

/((OOT);a;)^v = \A\'^sii?eB + |5p cos^ ± sin20B sin(3^ T 5)- (4.31) 

which is the same as Eq. (14. lip except for reversing the sign for the term proportional to 
sin(3\l/ =F 5)- The scattering intensity for G = (002) is given by Eq. (14.31 p with replacing 
sin(3^T(5) by sin(3^±5). 



2. Circular polarization 

The [/i, U2, and for G = (OOT) are given by the same forms as Eqs. (HTTD -l liTTp 
with reversing the signs of the terms proportional to cos(3\E' =F 5) and sin(3\l' =F 5). Thereby, 
the scattering intensity with Pi = is expressed as 



/((OOl); u) = /o((001); co) + /i((001); co) sin(3^ T S), 



(4.32) 



with 



/o((001);a;) 



20 ,1 .|2l + sin^^B 



\BY cos' en + \A 



X (1 + sin^ Ob ± 2P2 sin ^b + P3 cos^ ^b), 
/i((OOT);w) = |AP|cos0B(-i^2COs2^B ±2P3sin^B). 



(4.33) 
(4.34) 

Equation fl4.33p for Jo is the same as Eq. (7) in Ref. |5| (see also the errata), if |Ap and |Pp 
are identified to be proportional to and T^. Equation (I4.34p for Ji is similar to Eq. (8) 
in Ref. 5| and to Eq. (11) in Refs. jlll|l6l |. but different from them concerningchirality 
dependence. In addition, the possible phase shift 5 is absent in Eq. (6) in Ref. j5| and in 



Eq. (9) in Ref. (ll|. 

Finally, Ui, U2, and for G = (002) are given by the same forms as those for 
G = (001) with replacing 3\E' =F 5 by 3\E' ± 5 and with reversing the signs in the expression 
of U2 for G = (001). Hence, the scattering intensity with Pi = is expressed as 



/((002); Lo) = /o((002); to) + /i((002); to) sin(3^ ± 5), 



(4.35) 
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with 



/o((002);w) = l^l^cos^^B + 1^1 



,1 + sin2 Ob 



X (1 + sin^ Ob T 2P2 sin Ob + P3 cos^ ^b), 
Ji((002);w) = I 1 cos ^6(^^20052 ^B±2P3 sin ^b)- 



(4.36) 
(4.37) 



V. BOND-ORBITAL MODEL ANALYSIS 

We employ the bond-orbital model developed by Harrison^^ to describe the electronic 
structures of a-quartz and a-berlinite. We explain the model for a-quartz by following Ref. 
[1^ . The extension of the model to a-berlinite is straightforward by replacing one of the Si 
atoms by an Al atom and another one by a P atom in the Si-O-Si bond. 

Consider the Si-O-Si bond shown in Fig. HI There exist four such bonds for each Si atom. 
Let us introduce the Si sp^-orbitals, |/i(Si)). For the O atom, we construct the oxygen 
hybrids to the left and right, 



\h{0)i) = cosr/ls) + sinr7(— cos6'|p2) — sin^^lp^)), 
\h{0)r) = cosr^ls) + sin?7(cos6'|p2) — sin6'|pa;)). 



(5.1) 
(5.2) 



where \s) stands for the 2s state, and \px), \py), \pz) stand for the 2p states of O, respectively. 
Imposing the orthogonality condition {h{0)r\h{0)i) = 0, we have cos^i] = (1 — tan^^)/2 
with 6 = 18° being the bending angle of the bond. We may call the states by sp^'^^-hybrid. 
We also construct the lone-pair orbital, which we write as 



\hip) = cos r]ip\s) + sin r]ip\px). 



(5.3) 



Imposing the orthogonality conditions to \h{0)i) and \h{0)r), we have cos'^ rjip = tan^ ^. 
This state as well as the oxygen vr orbital, |7r(0)), are assumed to have no coupling to other 
states. With the bases, |/ii(Si)), \h{0)i), \h{0)r), |/i2(Si)), the Hamiltonian matrix for each 
bond may be represented as 

/ eu -V2 \ 

-V2 e,,(0) -V,- 
e,(0) -V2 

V -V2 eh2 ) 



bond 



(5.4) 
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\ip) 




FIG. 4: Si-O-Si bond for a-quartz. Si atoms are replaced by Al and P atoms for a-berlinite. 



e;,(Si) = -9.38 eV 
e/,(Al) = -6.96 eV 
e;,(P) = -11.96 eV 



TABLE I: Parameters for Si02 and AIPO4. 

e^(0) = -24.46 eV 
eip{0) = -18.55 eV 
e^(0) = -16.77 eV 
(9 = 18 deg 



^2(81-0) = 9.47eV 
y2(Al- O) = 8.20 eV 
V2{F-0) = 10.62 eV 
Vi- = 7.74 eV 



where eti = eh2 = e/i(Si) is the energy of the sp^ state of Si. The e/i(0), the hybrid energy 
of O, is given by ^h{0) = cos^ V^s{0) + sin^ rjep^O), with 6^(0) and 6^(0) being the energies 
of 2s and 2p states of O. The covalent energy V2 is given by 



1 /I a/S \ a/S 

V2 = -- COS r]Vssa + I - sin ?7 + — COS ?7 j V,p„ + — sin r]Vpp„, 



(5.5) 



where Vssa, Vspa, and Vppa are the Slater-Koster parameters.— The Vi- represents the cou- 
phng between the two oxygen hybrids given by 



Vi_ = (ep(0)-e,(0))cos2r7. 



(5.6) 



Table H] shows the parameter values for a-quartz and a-berlinite. Most of them are taken 



from Table 12-1 in Ref. 



12|. The eft(Al) and e/,(P) are taken from Table 1-1 in Ref. |l2|. 



The covalent energy V2 for Al-0 and P-0 are estimated from V2 for Si-0 by assuming the 
so-called dependence. The scattering intensities discussed later have been checked to be 
insensitive to the choice of parameter values. 



We diagonalize the Hamiltonian matrix to obtain the eigenvalues specified as e^^"* < e 



_(o) 



< 



.(0) 



< e 



(0) 

4 5 



and the corresponding eigenstates denoted by {i 



■4). In the ground 



15 



state, four electrons are occupied on the lowest two levels, \ipi ) and \iIj2 )) and furthermore 
four electrons are occupied on \hip) and |7r(0)). 
The El transition may be expressed as 

^^\9) = (5.7) 

where Xa is the dipole operator on the Si site. T = f R^p{r)rRis{r)r'^(\r with Rzpij") and 
-Ris(r) being the radial wave function of 3p and Is states. The aiso- is the annihilation 
operator of the Is electron with spin a, and ci\fj{tj is the creation operator of electron on 
orbital \iljf^)i in the £-th bond, where four bonds of Si-O-Si (or Al-O-P) are distinguished by 
£. The cos ■* denotes the directional cosine of the sp^ hybrid of the £-th bond with respect 
to the which is determined from the positions of oxygen atoms specified in Sec. [Tll 

The sum over i should be restricted within the unoccupied states, i.e. i = 3,4. Coefficient 
hf'^ is defined by \ipf'^)e = hf\hi)i + ■ ■ ■ , which is independent of £. 

In the intermediate state, an attractive potential from the Is-core hole is working, which 
is known to modify the absorption coefficient as a function of photon energy cu.— Assuming 
that the core-hole potential works within the Si or Al site, we change {hi\H\hi) from ehi 
to e/ji + Kore- We tentatively put V^orc = —6 eV. As shown below, the dependence on 
polarization, chirality, and scattering vector is unaltered by the presence of Kore for fixed cu, 
although the spectral shape is modified as a function of uj. We diagonalize the Hamiltonian 
matrix to obtain the eigenvalues e^^^ < 63^^ < e~^^ < 64^^ and the corresponding eigenstates 
iV^i^^) (^ = 1; ■ ■ '4). The intermediate state is constructed by distributing five electrons on 
these energy levels in one of the four bonds (other three bonds are occupied four electrons). 

The dominant contributions come from the intermediate states that four electrons occupy 
lowest two states and one electron occupies higher state. Thus, neglecting the so-called 
shake-up states in the intermediate state, we consider the overlap between the intermediate 
state and the ground state, and obtain the local dipole-dipole correlation function, 

[pt\e)U = \T'Q{e)Y,^o.e(^ cos<\ (5.8) 

where 

Q(e) = \OX'{\Oz,,hf + 0,,,hf\H{e - eeore - - e^'^) 

+ \0,^,hf + r5(e - e,ore - 6°, - ei'^)}, (5.9) 
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with 



Og = 


{#■!" 
(4" 






*4°'> 
4°') 






{•^•!" 






4°') if^'l^') 


On,n' 








4°') {4"l€') 






^■!°'> 




4°') {4"i4°') 


^ex = < 


Hi" ■ 


f4" 


- 


,(0)^ 


. 



(5.10) 

(5.11) 
(5.12) 

We notice from Eq. 05.81) that aU the components a(e), 6(e), c(e), d{e), are proportional to 
(5(e), and that their relative ratios are determined by the directions of Si-0 or Al-0 bonds. 
Hence we obtain A{u) and B{uj) from Eqs. (13.131) and (I3.14p . which are proportional to 



:de. 



u — e + iT 

Accordingly B*{uj)A{u) is real for both a-quartz and a-berlinite, indicating that the phase 
shift 6 is generally or vr, independent of u. These simple result would be modified, if the 
coupling between bonds is taken into account. 

The absorption coefficient C{uj) for photon energy u may be expressed as 

r/vr 



n\Xr 



(a; + e,-e„)2 + r2- 



(5.13) 



Therefore we obtain from Eqs. (13.61) and (15.81) 

r/vr 



C{uj) oc J Tr/3(j;e) 



(tu - e)2 + r2 



de = 2T^- 



Q(e) 



:de. 



[CO 



(5.14) 



VI. CALCULATED RESULTS IN COMPARISON WITH EXPERIMENTS 



A. a-quartz 

Figure [5] shows the absorption coefficient calculated from Eq. (I5.14p . in comparison with 
the experiment.- The first and second peaks arise from the transition from the Is state to 
the unoccupied state and to respectively. The core-hole potential makes the 

first peak intensity increase. In actuality, there exist more p-sjTiimetric states forming band 
states, which make the intensities spread above the K edge. The core-hole potential gives 
rise to the intensity transfer from the high energy region, resulting in sharpening the peak 
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FIG. 5: Absorption coefficient as a function of photon energy for Si-i^-edge in a-quartz. The sohd 
and broken Hnes are the calculated results with and without the core-hole potential. The origin of 
photon energy is set to the main peak position. Inset displays the experimental absorption data 
taken from Ref. ^. 

at the edge, as shown in Ref. jl8|. The states constituting the peak are close to the localized 
states constructed by the sp^-orbitals. Therefore the present model could describe rather 
well the states of the peak at the edge. Note that there exists a shoulder in the experimental 
curve of the absorption coefficient.- Such a peak is, however, not seen in the experimental 



and the theoretical curves in Ref. 



18| . The origin of this shoulder is not clear. 



Note that the RXS intensity as a function of u is proportional to 



uj — e + iV J yj oj — e + iV 

Although it is not shown here, this quantity has a large peak at uj giving the large peak in 
the absorption coefficient, in agreement with the experiment .- 

We concentrate our attention on the spectra at the uj giving the main absorption peak 
in the following. Figure [6] shows the RXS intensity as a function of azimuthal angle ^ for 
G = (001). Although G is defined by kj — kj in the experiment, which is opposite to ours 
(see the errata in Ref. {4]), we present G's in our definition. The Stokes parameters are set 
to be Pi = 0, P2 = ±0.95, P3 = —0.31 in accordance with the experiment.- For G = (001), 
sin6'B = 0.625, cos6'b = 0.781. The calculated intensities are larger for RCP (P2 = 0.95) 
than for LCP (P2 = —0.95) in No. 152, while the former is smaller than the latter in No. 154, 
consistent with the experimental curves shown in panels (b) and (d).- The intensities for 
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FIG. 6: RXS intensity from a-quartz as a function of azimuthal angle ^ for G = (001). The 
u) is fixed at the value giving the absorption peak. Pi = 0, P3 = —0.31, P2 = 0.95 (RCP), 
and P2 = —0.95 (LCP). Panels (b) and (d) show the experimental data in No. 152 and No. 154, 
respectively, which are taken from Ref. J]. 

LCP in No. 152 as well as for RCP in No. 154 are, however, too small in comparison with the 
experiment. We hope that the absorption correction to the experimental data, if it were not 
made yet, as well as a careful subtraction of the background might improve the discrepancy. 
As regards the oscillation terms, we notice from the general expressions (Eqs. f l4.32p - fl4.34p ) 
that they take the form of — ai sin(3\E' — 6) for RCP in No. 152, and ai sin(3\E' + 6) for 
LCP in No. 154 with ai a positive number, and that they take the form of 61 sin(3\[^ — 6) 
for LCP in No. 152, and — 61 sin(3\[^ + 6) for RCP in No. 154 with bi a positive number. 
The amplitude ai is much larger than the amplitude 61; their ratio is given by &i/ai = 
(IP2I cos^ 9b + 2P3 sin^^B)/(|-P2| cos^ 6'b — 2P3 sin6'B) ~ 0.2, which is independent of the model. 
The experimental curves seem to belong to these forms with the phase shift 5 = 27r/3 ~ tt. 
The present calculation within the bond-orbital model gives the phase shift 6 = it. 
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FIG. 7: RXS intensity from a-quartz of No. 154, for G = (001) and (001). The u is fixed at the 
value giving the absorption peak. Pi = 0, P3 = -0.31, P2 = 0.95 (RCP), and P2 = -0.95 (LCP). 
Lower panel shows the experimental curves reproduced from Ref. [4|. 

Figure [7] shows the RXS intensities for both G = (001) and G = (OOT) in No. 154. 
According to the general expressions (Eqs. fl4.18p and (14.321) ). the oscillation terms are 
proportional to sin(3\E' + 6) for both G = (001) and (001), regardless of RCP or LCP, in 
No. 154. This means that all curves have to be maximum or minimum at the same \E'- values. 
This requirement seems not to be satisfied in the experimental curves, where the maximum 
and minimum positions are somewhat different, indicating that the phase does not seems to 
has a definite value. 



B. a-berlinite 

Figure M shows the absorption coefficient calculated from Eq. (I5.14p , in comparison with 
the experiment.- Without taking account of the core- hole potential, we have the second- 
peak intensity larger than the first-peak one, which in fact is different from a-quartz. The 
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FIG. 8: Absorption coefficient as a function of plioton energy for Al-X-edge in a-berlinite. The 
solid and broken lines are the calculated results with and without the core-hole potential. The 
origin of photon energy is set to the main peak position. The inset displays the experimental 
absorption data taken from Ref. 5|. 

core-hole potential makes the first-peak intensity larger than the second one, in agreement 
with the experiment .- 

We concentrate our attention on the spectra at the uj giving the absorption peak. Figure 
|9]shows the RXS intensity for G = (001), where sin6'B = 0.361, cos = 0.932. The Stokes 
parameters are set to be Pi = 0, P2 = ±0.95, P3 = 0.30, in accordance with the experiment,- 
where the value of P3 here is opposite in sign to the case of a-quartz.- The average intensity 
for RCP is larger than for LCP in No. 152, while the former is smaller than the latter in 
No. 154, which dependence is the same as in a-quartz. The oscillation terms take the same 
forms as for a-quartz; the ratio of the amplitudes are given by 61/ai ~ 1.71. The present 
calculation gives the phase shift 5 = vr. The calculated curves are in good agreement with 
the experimental curves shown in panel (c).- We would like to emphasize that there exist 
no adjustable parameter in the present calculation. 

Figure HO] shows the RXS intensity for G = (002), where sin^B = 0.723, cos^b = 0.691. 
The Stokes parameters are the same as for G = (001). The intensity for LCP is larger 
than that for RCP in No. 152, while the former is smaller than the latter in No. 154. This 
dependence on polarization is opposite to that in G = (001). According to the general 
expressions (Eqs. fl4.35p and f l4.37p ). the oscillation terms take the form of 02 sin (3 + 5) 
for RCP in No. 152, and — a2 sin(3\[^ — 5) for LCP in No. 154 with 02 a positive number. 
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FIG. 9: RXS intensity from a-berlinite as a function of azimuthal angle ^, for G = (001). The 
uj is fixed at the value giving the absorption peak. Pi = 0, P3 = 0.30, P2 = 0.95 (RCP), and 
P2 = —0.95 (LCP). Panel (c) shows the experimental data for No. 154, taken from Ref. [5]. 

while they take the form of — 62sin(3\l' + 6) for LCP in No. 152, and 62 sin(3^ — 6) for 
RCP in No. 154 with 62 a positive number. The ratio of the amplitude is given by 62/02 = 
(IP2I cos^ 9b — 2P3 sin 6'b)/(|-P2| cos^ 9b + 2P3 sin^^e) ~ 0.02, corresponding to the nearly flat 
curves for LCP in No. 152 and for RCP in No. 154. The calculated curves are in good 
agreement with the experimental curves shown in Panel (c). 




VII. CONCLUDING REMARKS 



We have analyzed the RXS spectra at Si and Al K-edges on forbidden spots in chiral 
materials, a-quartz and a-berlinite. Summing up the local scattering matrices which satisfy 
the symmetry requirement, we have derived the general expression of scattering matrix on 
forbidden spots with El-El process. We have obtained the oscillation term as a function 
of azimuthal angle, which has the phase shift and the amplitude with chirality dependence 
different from the previous studies. We have evaluated the parameters undetermined by the 
symmetry argument alone on the basis of the underlying electronic structures given by the 
bond-orbital model. With such evaluation, the scattering matrix is completely determined. 
We have calculated the spectra depending on polarization, chirality, and scattering vector 
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FIG. 10: RXS intensity from a-berlinite as a function of ^, for G = (002). The u is fixed at the 
value giving the absorption peak. Pi = 0, P3 = 0.30, P2 = 0.95 (RCP), and P2 = -0.95 (LCP). 
Panel (c) shows the experimental data for No. 154, taken from Ref. j5l|. 

in agreement with the experiments, although some discrepancies remain in the average 
intensities and the phase shift of oscillation in a-quartz. The spectra for a-berlinite have 
reproduced particularly well the experiment. It should be emphasized that this result is 
obtained without adjustable parameters. 

We have obtained the phase shift of oscillation vr in both a-quartz and a-berlinite. For 
a-berlinite, the phase shift vr is consistent with the experiment. Note that the same phase 
shift is also observed in the chiral metal Te.— On the other hand, the situation of a-quartz 
is different, where the phase shift looks deviating from vr with depending on polarization 
and scattering vector (see the lower panel in Fig. [7j). The present formula based on the 
symmetry requirement does not allow such dependence but allows dependence on photon 
energy. Since both materials belong to the same covalent-bonding with similar electronic 
structures, such difference in the spectra between a-quartz and a-berlinite is puzzling to 
us. To clarify the origin of the phase shift and to obtain better agreement with experiment, 
it may be necessary to take account of the coupling between the bond orbitals, or more 
precisely, the band effects. It might be necessary to consider the absorption correction to 
the experimental curves, if it were not made, for more quantitative comparison. 

Finally we comment on the effect of the E1-E2 process. Since the space inversion symme- 
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try is broken around Si sites in a-quartz and around Al sites in a-berlinite, the p-symmetric 
states could mix with the (i-symmetric states, and therefore the second-order process using 
both the El transition and the electric-quadrupole (E2) transition could take place. This 
E1-E2 process is known to be important for x-ray absorption as well as RXS at the pre-K- 
edge in the transition-metal compounds such as a-Fe203,^° K2Cr04,^^ magnetite,— and 
GaFeOs-— These spectra have been analyzed based on the symmetry or by taking account 
of the microscopic electronic structures.— Now, for a-quartz, it has been proposed^'^'^i'ii 
that the phase shift is brought about by adding the scattering amplitude coming from the 
E1-E2 process to that from the El-El process. Since the former contribution is expected 
to be more than one order of magnitude smaller than the the latter, we think it unlikely to 
expect the substantial phase shift deviation from this mechanism. 
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